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Recently,  remarkable  advances  have  been  made  in  coupling  a  num¬ 
ber  of  high-Q  modes  of  nano-mechanical  systems  to  high-finesse 
optical  cavities,  with  the  goal  of  reaching  regimes  in  which  quan¬ 
tum  behavior  can  be  observed  and  leveraged  toward  new  applica¬ 
tions.  To  reach  this  regime,  the  coupling  between  these  systems 
and  their  thermal  environments  must  be  minimized.  Here  we  pro¬ 
pose  a  novel  approach  to  this  problem,  in  which  optically  levitating 
a  nano-mechanical  system  can  greatly  reduce  its  thermal  contact, 
while  simultaneously  eliminating  dissipation  arising  from  clamp¬ 
ing.  Through  the  long  coherence  times  allowed,  this  approach 
potentially  opens  the  door  to  ground-state  cooling  and  coherent 
manipulation  of  a  single  mesoscopic  mechanical  system  or  entan¬ 
glement  generation  between  spatially  separate  systems,  even  in 
room-temperature  environments.  As  an  example,  we  show  that 
these  goals  should  be  achievable  when  the  mechanical  mode 
consists  of  the  center-of-mass  motion  of  a  levitated  nanosphere. 

entanglement  |  optical  levitation  |  quantum  information 

One  of  the  most  intriguing  questions  associated  with  quantum 
theory  is  whether  effects  such  as  quantum  coherence  and 
entanglement  can  be  observed  at  mesoscopic  or  macroscopic 
scales.  As  a  first  step  toward  resolving  this  question,  recently 
much  effort  has  been  directed  toward  quantum  state  preparation 
of  high-Q  modes  of  nano-  and  micro-mechanical  oscillators — in 
particular,  cooling  such  modes  to  their  quantum  ground  state  (1). 
Reaching  a  regime  in  which  quantum  properties  such  as  entan¬ 
glement  (2)  emerge  is  not  only  of  fundamental  interest  but  could 
lead  to  new  applications  in  fields  such  as  ultrasensitive  detection 
(3,  4)  and  quantum  information  science  (5,  6).  To  reach  this  re¬ 
gime,  it  is  critical  that  the  thermalization  and  decoherence  rates 
of  these  systems  be  minimized  by  reducing  the  coupling  to  their 
thermal  reservoirs.  Thus  far,  this  has  necessitated  the  use  of  cryo¬ 
genic  operating  environments.  From  an  engineering  standpoint,  it 
would  also  be  desirable  to  reduce  the  dissipation  and  thermal¬ 
ization  rates  of  these  systems  through  their  clamping  and  material 
supports  (7),  so  that  these  rates  might  approach  their  fundamen¬ 
tal  material  limits  (8). 

Here  we  propose  a  unique  approach  toward  this  problem, 
wherein  the  material  supports  are  completely  eliminated  by 
optically  levitating  (9)  a  nano-mechanical  system  inside  a  Fabry- 
Perot  optical  cavity.  Indeed,  since  the  pioneering  work  of  Ashkin 
on  optical  trapping  of  dielectric  particles  (9)  (in  the  classical 
domain),  it  has  been  realized  that  levitation  under  good  vacuum 
conditions  can  lead  to  extremely  low  mechanical  damping  rates 
(10,  11).  We  show  that  such  an  approach  should  also  facilitate 
the  emergence  of  quantum  behavior  even  in  room-temperature 
environments,  when  the  particles  are  of  subwavelength  scale  such 
that  the  effects  of  recoil  heating  due  to  scattered  photons  become 
small.  As  a  specific  example,  we  show  that  the  center-of- 
mass  (CM)  motion  of  a  levitated  nanosphere  can  be  optically 
self-cooled  (12-14)  to  the  ground  state  starting  from  room 
temperature.  This  system  constitutes  an  extreme  example  of 
environmental  isolation  because  the  CM  motion  is  naturally 


decoupled  from  the  internal  degrees  of  freedom  in  addition 
to  being  mechanically  isolated  by  levitation.  In  this  case,  the 
decoherence  and  heating  rates  are  fundamentally  limited  by  the 
momentum  recoil  of  scattered  photons  and  can  be  reduced 
simply  by  using  smaller  spheres.  The  long  coherence  time  allowed 
by  small  spheres  enables  the  preparation  of  more  exotic  states 
through  coherent  quantum  evolution.  Here,  we  consider  in  detail 
two  examples.  First,  we  describe  a  technique  to  prepare  a 
squeezed  motional  state,  which  can  subsequently  be  mapped  onto 
light  leaving  the  cavity  using  quantum  state  transfer  protocols 
(15-18).  Under  realistic  conditions,  the  output  light  exhibits 
up  to  ~15  dB  of  squeezing  relative  to  vacuum  noise  levels,  poten¬ 
tially  making  this  system  a  viable  alternative  to  traditional  tech¬ 
niques  using  nonlinear  crystals  (19,  20).  Second,  we  show  that 
entanglement  originally  shared  between  two  modes  of  light  (21) 
can  be  efficiently  transferred  onto  the  motion  of  two  spheres  trap¬ 
ped  in  spatially  separate  cavities,  creating  well-known  Einstein- 
Podolsky-Rosen  (EPR)  correlations  (22)  between  the  mechanical 
systems.  Our  approach  of  optical  levitation  mirrors  many  suc¬ 
cessful  efforts  to  cool  (23,  24),  manipulate  (25)  and  entangle 
(26)  the  motion  of  atoms  and  ions  in  room-temperature  envir¬ 
onments.  At  the  same  time,  our  system  has  a  number  of  po¬ 
tential  advantages,  in  that  it  enables  direct  imaging  via  strong 
fluorescence,  exhibits  large  trap  depths,  and  has  a  relatively  large 
mass.  We  also  note  recent  related  experiments  involving  opto¬ 
mechanical  “fluids”  (with  a  continuous  excitation  spectrum  rather 
than  discrete  modes)  in  the  form  of  trapped,  ultracold  atomic  gases 
(27,  28). 

Beyond  the  examples  presented  here,  the  use  of  a  levitated 
device  as  an  opto-mechanical  system  could  provide  opportunities 
on  a  diverse  set  of  fronts.  For  instance,  it  should  allow  mechanical 
damping  to  approach  fundamental  material  limits,  potentially 
enabling  the  exploration  of  nanoscale  material  properties.  By 
levitating  systems  with  internal  vibrational  modes,  multiple 
modes  could  be  optically  addressed  and  cooled.  In  addition,  the 
CM  oscillation  frequency  can  be  tuned  through  the  trapping  in¬ 
tensity,  allowing  for  adiabatic  state  transfer  (29)  with  other  modes 
or  matching  spatially  separate  systems  for  optical  linking  and 
entanglement  generation  (30).  Furthermore,  this  paradigm  inte¬ 
grates  nano-mechanics  with  many  techniques  for  atomic  trapping 
and  manipulation,  which  can  be  further  extended  by  levitating 
systems  containing  an  internal  electronic  transition  (e.g.,  a  color 
center  within  a  nano-crystal)  (31).  Finally,  as  illustrated  by 
squeezed  light  generation,  engineering  mechanical  nonlinearities 
in  conjunction  with  quantum  state  transfer  yields  a  unique  means 
to  realize  nonlinear  optical  processes. 
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Optical  Forces  and  Noise  Acting  on  a  Dielectric  Sphere 

To  illustrate  our  idea,  we  consider  a  subwavelength  dielectric 
sphere  interacting  with  two  standing-wave  optical  modes  of  a 
Fabry-Perot  cavity  (Fig.  L4).  One  resonantly  driven  mode  pro¬ 
vides  an  optical  dipole  trap  for  the  sphere.  The  second  mode 
is  driven  by  a  weaker  “cooling”  beam,  assumed  to  have  a  nonzero 
intensity  gradient  at  the  trap  center,  which  provides  a  radiation 
pressure  cooling  force  (12-14).  We  discuss  the  cooling  mechan¬ 
ism  in  the  next  section,  whereas  here  we  focus  on  the  trapping 
potential  and  the  noise  forces  acting  on  the  sphere. 

The  trapping  beam  provides  a  gradient  force  similar  to  that 
used  to  “optically  tweeze”  small  dielectric  particles  (9).  Consid¬ 
ering  a  sphere  whose  radius  is  much  smaller  than  the  optical 
wavelength,  r  < c  2,  its  optical  response  is  like  that  of  a  point 
dipole  with  induced  dipole  moment  pind  =  aindE(x)  and  optical 
potential  Uopt(x)  =  -(l/4)(Reaind)|E(x)|2  (SI  Text).  Here  x  is 
the  CM  position  of  the  sphere  along  the  cavity  axis,  amd  = 
3e0V(^)  is  its  polarizability,  V  is  the  sphere  volume,  and  e  is  the 
electric  permittivity.  Taking  a  standing  wave  E(x)  =  Edcoskx 
{k  =  2 n/X),  to  lowest  order  near  an  antinode  the  potential  corre¬ 
sponds  to  a  harmonic  oscillator  with  mechanical  frequency 


«2'»Re'-lV'2 


pc 


e  +  2 


[1] 


where  I0  is  the  field  intensity  and  p  is  the  mass  density  of  the 
sphere.  The  total  trap  depth  is  U0  =  (3 10V /c)  Re  Typical  trap 
depths  and  oscillation  frequencies  are  plotted  in  Fig.  1C  and  D. 
For  all  numerical  examples,  we  take  material  properties  e  =  2, 
p  =  2  g/cm3  corresponding  to  fused  silica  and  an  operating  wa¬ 
velength  2=1  pm.  Frequencies  of  com/2n  ~  0.5  MHz  are  achiev¬ 
able  using  an  intracavity  intensity  of  /0~0.1W/pm2.  The 
imaginary  component  of  e  characterizes  optical  absorption,  which 
contributes  to  an  increase  ATint  in  the  internal  temperature  of  the 
sphere.  Assuming  a  value  corresponding  to  ~10  dB/km  propaga¬ 
tion  losses  in  bulk,  intensities  of  70  ~  10  W/pm2  can  be  sustained 
without  melting  the  sphere,  due  to  blackbody  reradiation  of  the 


absorbed  energy  (SI  Text).  We  believe  that  this  loss  value  is 
realistic,  given  that  even  lower  values  around  these  wavelengths 
have  been  observed  in  fused  silica  (32). 

The  dominant  noise  forces  acting  on  the  sphere  are  collisions 
with  a  background  gas  and  momentum  recoil  kicks  due  to  scat¬ 
tered  photons.  In  SI  Text ,  we  show  that  the  contributions  from 
shot  noise,  blackbody  radiation,  and  sphere  anisotropy  are 
negligible.  Furthermore,  the  CM  is  decoupled  from  the  internal 
degrees  of  freedom  and  the  sphere  effectively  has  no  internal 
structure  (as  opposed  to  molecules,  where  the  internal  configura¬ 
tion  can  affect  cooling  efficiency)  (33).  In  the  regime  in  which  the 
molecular  mean  free  path  exceeds  r,  the  background  gas  leads 
to  a  mean  damping  force  dp/dt  =  ~YgP/2  with  damping  rate 
yg/ 2  =  (8 /ji)(P/vrp),  where  P  and  v  are  the  background  gas 
pressure  and  mean  speed,  respectively  (34).  The  random  nature 
of  the  collisions  also  thermalizes  the  motional  energy  at  a  rate 
given  through  the  fluctuation-dissipation  theorem  by  dE/dt  = 
-yg(E  -kBT),  where  T  is  the  gas  temperature.  In  particular, 
the  characteristic  time  for  the  system  to  heat  by  one  phonon  start¬ 
ing  from  the  ground  state  is  rg  =  hcom/ygkBT.  Note  that  t~1  does 
not  necessarily  reflect  the  actual  collision  rate  between  the  sphere 
and  gas  molecules,  Rcoll  «  nPvr2 /kBT  (it  is  possible  for  a  single 
collision  to  be  quite  rare,  Rcoll  <^c  t~\  and  to  impart  several  pho¬ 
nons  at  once).  We  define  a  mechanical  quality  factor  Qg  =  com/yg 
due  to  the  background  gas  and  a  number  of  coherent  oscillations 
N% c  =  a>mTg /2*n  expected  before  the  energy  increases  by  a  single 
phonon.  For  a  sphere  of  radius  r  =  50  nm,  com/2ji  =  0.5  MHz, 
and  a  room-temperature  gas  with  P  =  10-10  Torr,  one  finds 
yg  ~  10-6  s-1,  Qg  ~  3  x  1012,  Nose  ^4x  104,  indicating  that  an 
ideal  nanosphere  can  be  essentially  decoupled  from  its  thermal 
environment. 

Photons  scattered  by  the  sphere  out  of  the  cavity  lead  to  heat¬ 
ing  and  decoherence  of  the  motion.  In  particular,  a  scattering 
event  entangles  the  mechanical  and  photonic  wave  functions  be¬ 
cause  the  phase  of  the  outgoing  photon  becomes  correlated  with 
the  position  of  the  scatterer.  Tracing  out  the  photonic  degrees  of 
freedom,  the  density  matrix  p  describing  the  motion  evolves  after 
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Figure  1.  A)  Illustration  of  dielectric  sphere  trapped  in  optical  cavity.  The  large  trapping  beam  intensity  provides  an  optical  potential  U0 pt(x)  that  traps  the 
sphere  near  an  antinode.  A  second  more  weakly  driven  cavity  mode  with  a  nonvanishing  intensity  gradient  at  the  trap  center  is  used  to  cool  the  motion  of  the 
sphere.  B)  Energy  level  diagram  of  mechanical  motion  (denoted  m)  and  cavity  cooling  mode  (ph).  The  mechanical  mode  has  frequency  comi  while  the  optical 
mode  has  frequency  oo2  and  linewidth  k.  Photon  recoil  induces  transitions  between  mechanical  states  | nm)  ->  \  (n  ±  1  )m)  at  a  rate  Rn^n± i  (^o^i  shown  by  dashed 
gray  arrow).  The  cooling  beam,  with  effective  optomechanical  driving  amplitude  Qm,  induces  anti-Stokes  scattering  that  cools  the  mechanical  motion  and 
allows  for  quantum  state  transfer  between  motion  and  light.  This  beam  is  also  responsible  for  weaker,  off-resonant  heating  via  Stokes  scattering.  O  Mechan¬ 
ical  frequency  com  as  a  function  of  trapping  beam  intensity.  For  all  numerical  results,  we  take  A  =  1  \im ,  p  =  2  g/cm3,  and  e  =  2 .  D)  Optical  trap  depth  U0  (in  K)  as 
functions  of  trapping  beam  intensity  and  sphere  radius  r. 
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a  scattering  event  as  p  ^  f dkP(k)elkrpe~lkr  (35).  Here  P( k)  is 
the  probability  distribution  of  the  scattered  photon  wavevector, 
and  r  is  the  CM  position  operator.  In  analogy  with  atoms  or  ions 
trapped  in  the  Lamb-Dicke  regime  (25),  when  the  particle  is 
trapped  on  a  scale  Ax  1  /k,  the  terms  e±lk  r  can  be  expanded 
to  first  order  in  k.  This  describes  transitions  between  consecutive 
harmonic  oscillator  levels  n  ->•  n  ±  1,  with  rates  Rn^n±i  = 
ysc(n  +  1/2  ±  1/2).  Considering  motion  only  along  the  x  direc¬ 
tion, 

Y SC  (2/5)(fi?r/ft?m)i?  sc  [2] 

where  eor  =  fik2/2pV  is  the  recoil  frequency  andRsc  is  the  photon 
scattering  rate.  A  result  identical  to  Eq.  2  holds  for  a  weakly 
excited,  trapped  atom  (36).  In  the  case  of  the  sphere,  Rsc  = 
48tt3  (fz^)2,  as  can  be  obtained  by  taking  the  power  radiated 

by  a  dipole  of  strength  pmd  and  dividing  by  the  energy  fleet  per 
photon.  We  emphasize  that  in  the  Lamb-Dicke  regime,  the  effect 
of  photon  scattering  on  motional  heating  and  decoherence  is 
completely  described  by  Eq.  2.  On  the  other  hand,  the  above 
transformation  of  the  density  matrix  p  is  completely  general. 
For  example,  when  the  motional  wave  packet  has  a  spatial  extent 
Ax  ~  2,  one  finds  then  that  a  single  scattering  event  can  destroy 
quantum  coherence  (37).  Analogous  effects  also  occur  due  to 
absorption  and  emission  of  blackbody  photons.  In  our  system, 
however,  such  effects  are  negligible  because  Ax  «;  l/k  (i.e., 
cor/com  <3C  1  in  Eq.  2),  and  the  emission  rate  of  blackbody  photons 
is  much  smaller  than  Rsc. 

It  is  convenient  to  define  a  dimensionless  parameter 

Ajr2  £  -  1 

tp  =  Ysja)m=——(v/xi)  [3] 

5  e  +  2 

to  characterize  the  strength  of  photon  recoil  heating,  which 
primarily  depends  on  the  sphere  volume  compared  to  the  cubic 
wavelength.  This  scaling  reflects  the  fact  that  the  scattered  power 
and  dipole  force  scale  like  p2nd  and  pind ,  respectively.  We  empha¬ 
size  that  the  effect  of  recoil  heating  can  be  reduced  by  simply 
using  smaller  spheres.  In  the  limit  that  optical  scattering  is  the 
dominant  heating  process,  the  expected  number  of  coherent  os¬ 
cillations  is  Aqsc  =  1  /( 2^)  oc  A3 /V.  We  will  find  that  cj)  naturally 
appears  to  characterize  the  fundamental  limits  to  observing  quan¬ 
tum  behavior  in  our  system.  For  a  sphere  of  radius  r  =  50  nm, 

e  =  2,  and  A  =  1  pm,  one  finds  that  N^c  ~  150.  Comparing 
with  background  gas  collisions  at  P  =  10-10  Torr  and  com/2n  = 
0.5  MHz,  recoil  heating  dominates  Nosc  for  sphere  sizes 
r  >  10  nm.  Reaching  the  regime  Aosc  »  1  implies  that  the  sphere 
can  coherently  evolve  for  many  oscillation  periods  after  any  cool¬ 
ing  mechanisms  are  turned  off,  which  makes  this  system  a  pro¬ 
mising  candidate  for  observing  coherent  quantum  phenomena. 

Finally,  we  remark  thatRsc  can  be  very  large  ( Rsc  ~  1014  s-1  for 
I o  =  1  W /pm2  and  r  =  50  nm)  compared  to  atoms  or  ions,  which 
enables  direct  imaging.  The  large  scattering  is  due  to  the  large 
intensities  and  the  linear  response  of  the  sphere  (it  does  not 
saturate  like  an  atom  or  ion),  as  opposed  to  the  system  behaving 
as  a  lossy  element  in  the  cavity.  The  contribution  to  the  cavity  loss 
rate  is  ksc  =  l2n;2co(V2 /A3Vc)(^=k)2,  where  Vc  is  the  cavity  mode 
volume,  and  ksc  is  typically  much  smaller  than  the  natural  cavity 
linewidth  k. 

Cooling  the  Center-of-Mass  Motion  to  the  Ground  State 

We  now  describe  the  optical  cooling  effect  of  the  weaker,  second 
cavity  mode  (denoted  mode  2).  For  concreteness,  we  assume 
that  the  sphere  is  trapped  near  the  antinode  x  =  0  of  cavity  mode 
Ex  ex  cos  kpc  and  that  the  second  mode  has  spatial  profile  E2  cx 
cos(k2x  -  #/4)  such  that  the  intensity  gradient  is  maximized.  The 
total  Hamiltonian  of  the  system  is  given  in  a  rotating  frame  by 


H  =  -h8xa\ai  -  hgi(cos2k1x  -  1  )a\ax  -  h82d\a2 

p2 

-  hg2  cos  2  (k2x  -  n/A)a\a2  +  ^ 

pjcy 

+  [(fl\  +  CL\)  +  +  ^2)]'  [4] 

Here  p  and  x  are  the  momentum  and  position  operators  of  the 
CM,  at  is  the  photon  annihilation  operator  of  cavity  mode  i, 
and  Q  and  £2y/2£7  are  the  driving  amplitudes  of  modes  1  and 
2,  respectively.  8t  is  the  detuning  between  the  driving  field  and 
mode  frequency  when  the  sphere  sits  atx  =  0.  The  opto-mechan- 
ical  coupling  strengths  gt  =  characterize  the  position- 

dependent  frequency  shifts  due  to  the  sphere  (SI  Text),  where 
Vci  and  coi  are  the  mode  volume  and  resonance  frequency  of 
mode  i.  To  simplify  notation,  we  assume  that  modes  1  and  2  have 
similar  properties,  cox  «  co2  =  00,  etc.  In  addition  to  the  evolution 
described  by  H,  the  system  also  exhibits  cavity  losses  and  the  me¬ 
chanical  noise  described  previously. 

Expanding  the  opto-mechanical  coupling  term  of  mode  2 
around  x  =  0,  hg  cos  2(kx  -  n/A)a\a2  «  2hgkxa\a2,  one  finds  a 
linear  coupling  in  the  sphere  position  analogous  to  the  effect 
of  radiation  pressure  on  a  moving  mirror  of  a  Fabry-Perot  cavity 
(13).  Physically,  the  motion  induces  changes  in  the  detuning  and 
thus  the  intracavity  field  amplitude,  while  the  lag  in  the  cavity 
response  enables  the  field  to  do  work  (cooling)  on  the  sphere. 
Following  the  techniques  of  ref.  13,  to  calculate  the  cooling  rate 
we  first  apply  shifts  to  the  operators,  at  ->•  at  +  at  andx  ->  x  +x0, 
where  at  andx0  «  £/k  [£  «  k2C / (k2  +  4 8\)\  are  mean  values  of  the 
quantum  fields.  Here  we  have  defined  2f  =  \a2/ax\2  as  the  ratio 
of  intracavity  intensities  of  modes  1  and  2  and  assumed  that 
mode  1  is  driven  on  resonance  (^  =  0).  To  lowest  order  in  f ,  field 
mode  1  (mode  2)  is  purely  responsible  for  trapping  (cooling). 
Subsequently  tracing  out  the  cavity  degrees  of  freedom  yields 
equations  for  the  mechanical  motion  alone.  In  particular,  to 
lowest  order  in  f  and  for  82  <  0,  the  cooling  laser  provides 
a  net  cooling  rate  T  =  _Ropt>_  -  RovU+  =  kQ.1, [((<52  +  com)2+ 
(k/2)2)-1  -  ((<52  -  oomY  +  (k/2)2)-1]  (SI  Text),  where  i?opW  de- 
note  the  anti-Stokes  (cooling)  and  Stokes  (heating)  scattering 
rates  (see  Fig.  IB).  Here  klm  =  2gkxm\ai\y/2Z,  is  the  effective 
opto-mechanical  driving  amplitude  (see  Fig.  IB)  and  xm= 
^h/2mcom.  Validity  of  these  perturbative  results  requires  £lm 
<k,  com  and  £  <  1 . 

In  the  realistic  limit  that  background  gas  collisions  are  negli¬ 
gible,  the  steady-state  phonon  number  is  («y)  «  hf  +  ys C/F,  where 
hf  =  Ropt,+/ F  is  the  fundamental  limit  of  laser  cooling  (13).  It  is 
minimized  when  82  =  -(1/2)  sJk2  +  Aco2m.  In  particular,  when 
sideband  resolution  is  achieved  (com  >  k),  hf  m in  «  ( K/Acom )2  <^c  1, 
indicating  that  ground-state  cooling  is  possible  provided  other 
heating  mechanisms  are  made  sufficiently  small.  Considering 
the  limit  &>m  »  k  and  taking  the  maximum  cooling  rate  T  ~  k  con¬ 
sistent  with  the  perturbative  calculations,  using  Eq.  3  one  can 
then  rewrite  (nf)  as 

[5] 

The  last  term  on  the  right  corresponds  to  photon  recoil  heating. 
Eq.  5  is  minimized  for  K/com  =  2</>1/3,  in  which  case  (rif}min  = 
3(P2P/4  cx  ( r/A )2  <3C  1.  Thus,  one  sees  that  ground-state  cooling 
is  readily  attainable  (provided  that  f  <  1  can  be  simultaneously 
satisfied).  Physically,  the  optimum  value  of  K/com  balances  good 
sideband  resolution  and  excessive  heating  when  large  intensities 
are  used  to  increase  com. 

To  illustrate  these  results,  we  consider  a  fused  silica  sphere  of 
radius  r  =  50  nm  and  com/(2ji)  =  0.5  MHz  levitated  inside 
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a  cavity  of  length  L  =  1  cm  and  mode  waist  w  =  25  pm  (Vc  = 
( n/4)Lw 2).  In  Fig.  24  we  plot  the  minimum  obtainable  («y) 
(black  curve)  as  a  function  of  cavity  finesse  &  =  jic/kL ,  assuming 
negligible  gas  collisions  and  subject  to  the  constraints  2f,  £2m//c, 
^m/^m  <  1/2  and  optimized  over  detuning  <52.  For  low  cavity  fi¬ 
nesse  the  cooling  is  nearly  limited  by  sideband  resolution  (hf  m in, 
red  curve),  and  the  ground-state  regime  («y)  <  1  can  be  reached 
with  a  minimum  finesse  of  &  ~  104.  A  minimum  of  (rif)  ~  0.01  is 
reached  at  a  finesse  of  &  ~  3  x  105,  with  a  maximum  cooling  rate 
ofr~/c~3xl05  s_1  (see  Table  1  for  a  summary  of  parameters). 
This  corresponds  to  a  final  temperature  of  Tf  ~  5  pK,  or  a  re¬ 
markable  compression  factor  ofr/7y~6xl07  relative  to  room 
temperature  T.  We  note  that  values  of  &  ~  106  have  been 
achieved  in  Fabry-Perot  cavities  at  similar  wavelengths  (38). 

Thus  far,  we  have  only  considered  the  motion  along  the  cavity 
axis.  With  additional  cavities  oriented  along  different  axes,  clearly 
it  is  possible  to  achieve  3D  ground-state  cooling  of  the  CM  mo¬ 
tion.  Moreover,  due  to  the  high  rate  of  Rayleigh  scattering,  feed¬ 
back  cooling  can  be  readily  applied  for  transverse  localization. 
However,  as  shown  in  SI  Text ,  the  transverse  motion  need  not 
be  reduced  below  the  level  set  by  the  ambient  environment.  Spe¬ 
cifically,  one  only  requires  that  the  transverse  position  uncer¬ 
tainty  Ay  be  small  compared  to  the  beam  waist  w,  as  the 
transverse  motion  introduces  errors  of  order  (A y/w)2  (e.g.,  in 
the  fractional  decrease  in  cooling  efficiency).  For  the  parameters 
of  Table  1,  thermal  motion  at  T  =  300  K  yields  only  (A y/w)2  ~ 
10-2  due  to  the  large  trap  depth,  while  recoil  heating  leads  to  a 
position  uncertainty  A y/w  &  y/(p(Drt  in  the  absence  of  any  trans¬ 
verse  cooling  mechanism.  Note  that  result  of  recoil  heating  is  in¬ 
tensity-independent,  since  an  increase  in  intensity  leads  to  equal 
increases  in  the  photon  scattering  rate  and  the  transverse  confine¬ 
ment  of  the  optical  potential.  One  finds  that  recoil  heating  is  a 
very  weak  effect,  as  a  time  /~5,000  s  is  required  to  reach  the 
regime  (Ay/w)2  ~  10-2. 

Motional  Entanglement  and  Squeezed  Light  Generation 
Using  Quantum  State  Transfer 

A  number  of  related  schemes  have  been  proposed  for  quantum 
state  transfer  between  light  and  the  motion  of  atoms  (15,  16)  or 
nano-mechanical  systems  (17,  18).  In  our  system,  the  small  me¬ 
chanical  noise  and  ease  of  achieving  good  sideband  resolution  in 
principle  allow  state  transfer  to  be  realized  with  almost  perfect 
efficiency.  This  might  enable  light  with  nonclassical  properties 
to  be  mapped  onto  mechanical  motion  (39),  and  as  an  example, 
we  show  that  this  can  be  used  to  generate  EPR  correlations  be¬ 
tween  two  spatially  separate  spheres.  Moreover,  a  complemen¬ 
tary  process  can  be  realized,  where  a  nontrivial  mechanical 
state  (a  squeezed  state)  is  prepared  through  coherent  manipula¬ 
tion  and  subsequently  transferred  to  light  leaving  the  cavity.  The 
latter  case  illustrates  how  opto-mechanics  can  yield  a  novel 
nonlinear  optical  system. 


First  we  give  a  simplified  picture  of  quantum  state  transfer 
using  a  one-sided,  ideal  cavity  (where  all  losses  are  via  transmis¬ 
sion  through  one  cavity  mirror)  (40).  Specifically,  we  consider  the 
Heisenberg  equations  of  motion  in  a  rotating  frame  for  the  cavity 
cooling  mode  and  the  motion  (after  applying  the  shifts  described 
in  the  previous  section),  when  the  cooling  mode  is  driven 
resonantly  on  the  red  motional  sideband  (<52  =  -&>m), 

“  \ «2  “  ft  +  bfe2ia”‘t )  +  sjm 2M , 

^  =  (i/h)[He,b\  -  iQm(a2  +  a\e2l<°mt)  +  iF(t)e“°m‘.  [6] 

The  Hamiltonian  He  describes  any  external  forces  or  couplings 
applied  to  the  sphere  beyond  those  in  Eq.  4,  b  is  the  annihilation 
operator  corresponding  to  a  harmonic  oscillator  of  mass  m  and 
frequency  com,  and  <22  in  is  the  cavity  input  operator  associated 
with  losses.  F(t)  is  the  (Hermitian)  noise  due  to  photon  recoil, 
which  has  correlations  ( F(t)F{t '))  =  (f)0)m8{t  -t'),  and  we  assume 
all  other  noise  is  negligible.  Since  the  cavity  trapping  mode  ( a x) 
effectively  provides  a  harmonic  potential  and  can  otherwise  be 
ignored,  for  simplicity  we  will  omit  the  subscript  2  as  we  refer 
to  the  cooling  mode  in  future  discussions.  Temporarily  assuming 
that  the  nonsecular  terms  (e2l(°mt)  can  be  ignored  and  that 
the  mechanical  motion  evolves  slowly  on  time  scales  compared 
to  1/k,  one  can  adiabatically  eliminate  the  cavity  mode  to 
yield  a  «  -2 i(£lm/K)b  +  (2 /y/icja^,  and  db/dt «  (i/h)[He,b\- 
(r/2)S  —  iVrdin  +  iF(t)el(°mt,  where  T  =  4 is  the  cavity-in¬ 
duced  cooling  rate  in  the  weak-driving  limit  (Qm  <  k).  The  cavity 
output  is  related  to  the  input  and  intracavity  fields  through 
aOVLt  =  \[kci  -  ain  (40),  or  aout  «  —iy/Tb  +  ain,  which  states  that  the 
mechanical  motion  is  mapped  onto  the  outgoing  light.  Physically, 
the  cooling  process  converts  phononic  excitations  into  photonic 
excitations  that  leak  out  of  the  cavity.  Generally,  two  mechanisms 
will  degrade  state  transfer.  First,  F  adds  extra  noise  to  the  ideal 
state  that  one  is  attempting  to  transfer,  with  a  strength  character¬ 
ized  by  the  small  parameter  </>.  Second,  the  nonsecular  terms 
contribute  to  Stokes  scattering,  destroying  the  perfect  phonon- 
photon  correspondence,  with  a  strength  that  is  expected  to  be 
proportional  to  ( K/oom )2.  Given  that  </>,  (/ K/com )2  can  be  made 
small,  nearly  perfect  state  transfer  is  possible  in  principle.  We 
illustrate  this  with  two  examples,  entanglement  transfer  and 
squeezed  light  generation. 

Entanglement  Transfer.  Here  we  describe  how  EPR  correlations 
shared  between  two  modes  of  light  (21)  can  be  mapped  to  the 
motion  of  two  spheres  trapped  in  spatially  separate  cavities. 
Specifically,  we  define  quadrature  operators  for  the  input  light 

for  each  of  the  two  systems  (denoted  A ,  B ),  given  by  = 


Figure  2.  A)  Mean  phonon  number  (nf)  (black  curve)  versus  cavity  finesse  &  &  =  nc/xL)  under  optimized  cooling  conditions.  The  system  parameters  are  given 
in  Table  1 .  The  red  curve  denotes  hf  m in,  the  fundamental  limit  of  cooling  imposed  by  sideband  resolution.  B )  Solid  blue  curve:  optimized  EPR  variance  between 
two  levitated  spheres,  as  a  function  of  squeezing  parameter  e~2R.  System  parameters  are  identical  to  a).  Dashed  curve:  EPR  variance  in  limit  of  perfect  state 
transfer,  AEPR  =  e~2R.  Green  curve:  cavity  finesse  corresponding  to  optimal  EPR  variance.  Q  Optimized  variance  (AX^out)min  (in  dB)  of  squeezed  output  light 
from  an  ideal  cavity,  as  a  function  of  sphere  size. 
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Table  1.  Example  cooling  parameters  for  a  fused  silica  sphere  of  radius  r  =  50  nm  at  A  =  1 


Total  cavity  decay  Mechanical  Intracavity  intensities:  Photon  scattering 


Cavity 

length 

Beam 

waist 

Cavity 
finesse  & 

k/2ji,  scattering 
contribution  ksc/2ji 

frequency 

com/2n 

trapping,  cooling 
beams 

Internal  heating 

A7int 

rate  Rsc  (both 
beams) 

Phonon 
number  (nf) 

1  cm 

25  |im 

3  x  105 

5  x  1 04,  100  Hz 

0.5  MHz 

0.1,  0.05  W/nm2 

80  K 

1.6x1013  s"' 

0.01 

=  (flj^  -  a^)/i  for  j  =A,B.  A  similar  set 
of  operators  X±m,  X±out  can  be  defined  for  the  motion  and 

output  light,  by  replacing  a®  ->  ffl\  a^,  respectively.  Of  parti¬ 
cular  interest  is  the  case  where  the  two  input  fields  exhibit 
broadband  EPR  correlations  between  them, 

((*£L H  ±X%(a,)Y)/2  =  <  1.  [7] 

When  the  variances  satisfy  e~m  <  1,  the  two  modes  exhibit  cor¬ 
relations  below  vacuum  level  and  are  entangled  (41)  (for  concre¬ 
teness,  we  assume  the  other  combinations  of  quadratures  satisfy 

((X^(co)+Z^?n(ffl))2)/2  =  e2R).  Such  EPR  correlations  have 
been  observed  with  light  and  in  the  internal  degrees  of  freedom 
of  atomic  ensembles  (42),  but  have  yet  to  be  demonstrated  using 
mechanical  systems. 

To  proceed,  we  solve  Eq.  6  in  the  Fourier  domain  (including 
the  nonsecular  terms)  for  each  of  the  systems  for  the  correlations 
given  in  Eq.  7  and  He  =  0.  Generally,  the  nonsecular  terms  yield 
an  infinite  set  of  algebraic  equations  (coupling  frequencies  com  + 
2 ncom  for  integer  n),  which  given  com  ^klm,K  can  be  truncated  to 
good  approximation  at  n  >  1.  For  simplicity  of  analysis,  we 
assume  the  two  systems  have  identical  properties,  and  that  the 
cooling  rate  T  =  k.  However,  we  expect  our  results  should  qual¬ 
itatively  hold  provided  that  only  r,  com  of  the  two  systems  are 
properly  tuned,  which  can  be  easily  accomplished  by  adjusting 
the  trapping  and  cooling  beam  intensities.  One  can  then  show 
that  state  transfer  yields  the  following  joint  variances  in  the 
motion  (SI  Text), 

Aepr  =  ((X{fm(t)+X^m(t)Y)/2 

=  e~2R +X^-{3e2R +  2sinh2R)+^!E,  [8] 

16  OOfn  K 

As  expected,  Stokes  scattering  and  recoil  heating  contribute  to 
the  variance  by  amounts  ( K/com )2  and  (p(om/K ,  respectively.  This 
can  be  minimized  with  respect  to  K/com,  yielding  AEPR  min  =  e~m+ 
3(cf)/2)2^(3e2R  +  2sinh2R)1/3.  To  illustrate  these  results  we  plot 
AEpr, min  in  Fig-  27?  as  a  function  of  e~m,  taking  the  same  para¬ 
meters  as  in  Fig.  24.  For  the  moderate  values  of  e~2R  typically 
obtained  in  experiments  (21),  EPR  correlations  in  the  motion 
can  be  achieved  with  reasonable  cavity  finesse  &  ~  105. 

Squeezed  Light  Generation.  First  we  describe  a  technique  to  create 
a  mechanical  squeezed  state  and  then  derive  the  properties  of  the 
outgoing  light  upon  quantum  state  transfer.  Mechanical  squeez¬ 
ing  is  accomplished  by  adding  a  sinusoidally  varying  component 
to  the  intensity  of  the  trapping  beam,  which  yields  the  Hamilto¬ 
nian  of  a  parametric  amplifier  (43),  He  =  CnCo^x2  sin  2 eomt.  Here 
em  is  a  small  parameter  characterizing  the  strength  of  the  mod¬ 
ulation  of  the  trap  frequency.  As  one  approaches  the  threshold 
for  parametric  oscillation  (emcom  -»  T),  the  variance  in  one  quad¬ 
rature  of  motion  is  reduced  by  up  to  a  factor  of  2  (43). 

We  now  investigate  the  properties  of  the  outgoing  light  over  a 
narrow  frequency  range  near  the  cavity  resonance,  specifically 
considering  X±  out(co  =  0).  We  apply  similar  methods  as  above 
to  solve  Eq.  6  in  the  Fourier  domain.  Taking  the  limit  as  one 


approaches  threshold  and  T  =  k,  the  variance  in  the  output  light 
is  given  by  (SI  Text) 


A^i.out  (®  =  0) 


2  <t>a>m  5  k2 

K  1  6  (l)2n 


[9] 


Again,  an  optimum  value  of  K/com  oc  ft 1/3  maximizes  the  squeez¬ 
ing,  with  a  minimum  variance  of  (AA2  ,out)min  «2.04^2/3.  Thus 
for  small  sphere  sizes  the  noise  level  can  be  reduced  far  below 
that  of  the  vacuum  state  (corresponding  to  AA2  out  =  1).  A  plot 
of  (AA2  out)min  as  a  function  of  sphere  size  is  shown  in  Fig.  2C. 
For  r  =  10  nm  size  spheres,  one  finds  that  ~30  dB  of  noise  reduc¬ 
tion  relative  to  the  vacuum  state  can  be  obtained  using  an  ideal 
cavity  (note  that  for  a  background  gas  pressure  of  P  ~  10“ 10  Torr, 
additional  noise  arising  from  gas  collisions  is  negligible  down 
to  r  ~  10  nm). 

In  practice,  a  cavity  has  additional  scattering  and  absorption 
losses  that  limit  the  squeezing.  Starting  from  Eq.  9,  which  gives 
the  amount  of  squeezing  at  threshold  using  an  ideal  cavity  (with 
T  =  k),  we  model  cavity  losses  via  a  beam  splitter  transformation 
with  the  ideal  squeezed  light  and  vacuum  as  the  two  inputs.  The 
output  light  exhibits  reduced  squeezing  due  to  mixing  with  the 
vacuum,  given  by 


(A*iout(«  =  0))min=  (l  A) 


f  2(p<j)m  !  5  k2  '\ 
V  K  16a>m) 


[10] 


where  k',  k  denote  the  scattering/absorption  loss  in  the  cavity 
and  the  total  cavity  linewidth,  respectively.  In  the  relevant  re¬ 
gime  where  k'  /k  1,  we  can  approximate  1  -  k'  /k  «  1  and  the 
squeezing  is  optimized  for  the  choice  K/com  =  2(2/5) 1/3  +  k' / 

(2 (Bm))1/3,  for  which  (AW2  out)min  «  2.04(0  +  k1  / (2a>m))2'3 .  We 
now  must  choose  a  set  of  realistic  cavity  parameters  where  this 
optimized  squeezing  can  be  realized,  and  where  T  =  k  is  consis¬ 
tent  with  f  being  small.  As  an  example,  taking  a  cavity  length  and 
waist  of  L  ~  2  cm  and  w  ~  10  pm,  k'  corresponding  to  1  ppm 
losses  per  round  trip  (38),  and  sphere  parameters  r  =  50  nm  and 
com/(2n)  =  0.5  MHz,  we  find  that  T  =  k  corresponds  to  a  value 
f  ~  1/4,  which  yields  squeezing  of  ~15  dB  in  the  output  light. 

In  principle,  similar  techniques  also  apply  to  trapped  atoms  or 
ions.  However,  one  benefits  from  the  relatively  large  mass  m  of 
the  sphere.  Specifically,  approaching  threshold,  one  quadrature 
of  motion  becomes  infinitely  unsqueezed,  producing  a  large  posi¬ 
tion  uncertainty  Ax  (43).  At  the  same  time,  faithful  quantum  state 
transfer  requires  a  linear  opto-mechanical  coupling,  which  trans¬ 
lates  into  a  requirement  that  the  spatial  extent  Ax  of  the  motional 
state  of  the  nanosphere  be  well-localized  with  respect  to  the 
wavelegnth  A0  of  the  trapping  laser  and  the  wavelengths  As  for 
any  incoherent  scattering  events.  Taking  A0  =  As  =  2n/k  as  rele¬ 


vant  to  the  nanosphere,  we  require  that  kAx « k 


where  n  characterizes  the  average  phonon  number.  In  SI  Text, 
we  show  that  kAx  <  10-2  can  be  satisfied  with  a  sphere  even 
in  the  regime  of  ~30  dB  squeezing.  To  compare,  a  “typical”  atom 
of  mass  ma  =  100  amu  trapped  with  frequency  com/(2n) 
0.5  MHz  has  a  ground-state  wavepacket  of  size  xa  «  10  nm,  so 
that  Ax&xanx/2  «  320  nm.  Hence,  a  linear  coupling  cannot  be 
achieved  for  optically  trapped  atoms  for  the  requisite  n,  while 
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linear  coupling  is  possible  for  ion  traps.  However,  in  both  cases,  a 
single  scattered  photon  would  substantially  decohere  the  mo¬ 
tional  quantum  state  of  the  atom,  as  Ax  is  comparable  to  As. 

The  amount  of  squeezing  predicted  in  this  opto-mechanical 
system  compares  favorably  with  traditional  techniques  involving 
a  nonlinear  optical  crystal  pumped  inside  a  cavity,  where  squeez¬ 
ing  of -10  dB  has  recently  been  observed  (19,  20).  An  important 
limitation  in  these  experiments  appears  to  be  nonlinear  absorp¬ 
tion  of  the  crystal  (19),  which  currently  prevents  significantly 
greater  squeezing  even  if  all  other  imperfections  are  eliminated. 
On  the  other  hand,  the  mechanical  noise  in  our  system  in  prin¬ 
ciple  enables  squeezing  by  up  to  -30  dB,  limited  only  by  cavity 
losses  and  extraction  efficiency.  With  state-of-the-art  Fabry- 
Perot  cavities,  we  anticipate  that  up  to  -15  dB  of  squeezing  in  the 
output  light  can  be  achieved.  Even  higher  levels  should  be  pos¬ 
sible  by  coupling  the  levitated  particle  to  microsphere  cavities  to 
take  advantage  of  their  small  mode  volumes  and  remarkable 
quality  factors  (44),  provided  that  trapping  techniques  near  such 
structures  can  be  developed. 

Outlook 

An  optically  levitated  opto-mechanical  system  can  have  remark¬ 
ably  long  coherence  times,  which  potentially  enables  quantum 
phenomena  such  as  entanglement  to  be  observed  even  in  room- 
temperature  environments.  Combining  previously  demonstrated 
techniques  to  controllably  grow  small  particles  (45)  and  load 
and  manipulate  them  in  vacuum  (9,  46)  should  put  this  system 

1 .  Cleland  A  (2009)  Optomechanics:  Photons  refrigerating  phonons.  Nat Phys,  5:458-460. 

2.  Genes  C,  Vitali  D,  Tombesi  P  (2008)  Simultaneous  cooling  and  entanglement  of 
mechanical  modes  of  a  micromirror  in  an  optical  cavity.  New  J  Phys,  10:095009. 

3.  Mamin  HJ,  Rugar  D  (2001)  Sub-attonewton  force  detection  at  millikelvin  tempera¬ 
tures.  Appl  Phys  Lett,  79:3358-3360. 

4.  Maiwald  R,  et  al.  (2009)  Stylus  ion  trap  for  enhanced  access  and  sensing.  Nat  Phys, 
5:551-554. 

5.  Cleland  AN,  Geller  MR  (2004)  Superconducting  qubit  storage  and  entanglement  with 
nanomechanical  resonators.  Phys  Rev  Lett,  93:070501. 

6.  Rabl  P,  et  al.  (2009)  A  quantum  spin  transducer  based  on  nano  electro-mechancial 
resonator  arrays.  arXiv  0908.0316. 

7.  Hao  Z,  Erbil  A,  Ayazi  F  (2003)  An  analytical  model  for  support  loss  in  micromachined 
beam  resonators  with  in-plane  flexural  vibrations.  Sensor  Actuat  A-Phys,  109:156-164. 

8.  Wilson-Rae  I  (2008)  Intrinsic  dissipation  in  nanomechanical  resonators  due  to  phonon 
tunneling.  Phys  Rev  B,  77:245418. 

9.  Ashkin  A  (2007)  Optical  Trapping  and  Manipulation  of  Neutral  Particles  Using  Lasers:  A 
Reprint  Volume  with  Commentaires  (World  Scientific,  Singapore). 

10.  Ashkin  A,  Dziedzic  JM  (1976)  Optical  levitation  in  high  vacuum.  Appl  Phys  Lett, 
28:333-335. 

11.  Libbrecht  KG,  Black  ED  (2004)  Toward  quantum-limited  position  measurements  using 
optically  levitated  microspheres.  Phys  Lett  A,  321:99-102. 

12.  Braginsky  VB,  Vyatchanin  SP  (2002)  Low  quantum  noise  tranquilizer  for  Fabry-Perot 
interferometer.  Phys  Lett  A,  293:228-234. 

13.  Wilson-Rae  I,  Nooshi  N,  Zwerger  W,  Kippenberg  TJ  (2007)  Theory  of  ground  state  cool¬ 
ing  of  a  mechanical  oscillator  using  dynamical  backaction.  Phys  Rev  Lett,  99:093901. 

14.  Marquardt  F,  Chen  JP,  Clerk  AA,  Girvin  SM  (2007)  Quantum  Theory  of  Cavity-Assisted 
Sideband  Cooling  of  Mechanical  Motion.  Phys  Rev  Lett,  99:093902. 

15.  Zeng  H,  Lin  F  (1994)  Quantum  conversion  between  the  cavity  fields  and  the  center- 
of-mass  motion  of  ions  in  a  quantized  trap.  Phys  Rev  A,  50:R3589-R3592. 

16.  Parkins  AS,  Kimble  HJ  (1999)  Quantum  state  transfer  between  motion  and  light.  J  Opt 
B-Quantum  S  O,  1 :496-504. 

17.  Zhang  J,  Peng  K,  Braunstein  SL  (2003)  Quantum-state  transfer  from  light  to  macro¬ 
scopic  oscillators.  Phys  Rev  A,  68:013808. 

18.  Jahne  K,  et  al.  (2009)  Cavity-assisted  squeezing  of  a  mechanical  oscillator.  Phys  Rev  A, 
79:063819. 

1 9.  Takeno  Y,  Yukawa  M,  Yonezawa  H,  Furusawa  A  (2007)  Observation  of-9  dB  quadrature 
squeezing  with  improvement  of  phase  stability  in  homodyne  measurement.  Opt 
Express,  15:4321-4327. 

20.  Vahlbruch  H,  et  al.  (2008)  Observation  of  squeezed  light  with  10-dB  quantum-noise 
reduction.  Phys  Rev  Lett,  100:033602. 

21.  Yonezawa  H,  Braunstein  SL,  Furusawa  A  (2007)  Experimental  demonstration  of 
quantum  teleportation  of  broadband  squeezing.  Phys  Rev  Lett,  99:110503. 

22.  Einstein  A,  Podolsky  B,  Rosen  N  (1935)  Can  quantum-mechanical  description  of 
physical  reality  be  considered  complete?.  Phys  Rev,  47:777-780. 

23.  McKeever  J,  et  al.  (2003)  State-insensitive  cooling  and  trapping  of  single  atoms  in  an 
optical  cavity.  Phys  Rev  Lett,  90:133602. 

24.  Maunz  P,  et  al.  (2004)  Cavity  cooling  of  a  single  atom.  Nature,  428:50-52. 


within  immediate  experimental  reach.  Extending  the  ideas 
presented  here  should  open  up  several  other  interesting  possibi¬ 
lities.  First,  we  anticipate  that  our  techniques  to  achieve  environ¬ 
mental  isolation  and  use  optical  forces  to  engineer  the  mechanical 
motion  are  applicable  to  a  large  class  of  opto-mechanical  systems 
beyond  dielectric  nano-particles.  For  example,  by  introducing 
specifically  tailored  optical  potentials,  it  should  be  possible  to  pro¬ 
duce  nontrivial  dynamics,  such  as  nonlinear  motion.  In  addition, 
several  nano-particles  or  more  complex  nano-mechanical  systems 
with  internal  modes  could  be  levitated  and  coupled  together, 
for  the  purpose  of  entangling  multiple  degrees  of  freedom.  Finally, 
by  levitating  charged  or  magnetic  systems,  one  could  potentially 
realize  systems  analogous  to  ion  traps  (47)  or  facilitate  novel 
quantum  hybrid  architectures  (6). 

Note  added:  We  have  become  aware  of  a  recent,  similar 
proposal  to  optically  levitate  and  manipulate  a  nano-mechanical 
system  by  O.  Romero-Isart  et  al.,  in  arXiv:0909.1469. 
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